Giant vortices in combined harmonic and quartic traps 
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We consider a rotating Bose-Einstein condensate confined in combined harmonic and quartic traps, following 
recent experiments [V. Bretin, S. Stock, Y. Seurin and J. Dalibard, cond-mat/0307464]. We investigate numer- 
ically the behavior of the wave function which solves the three-dimensional Gross Pitaevskii equation. When 
the harmonic part of the potential is dominant, as the angular velocities f2 increases, the vortex lattice evolves 
into a giant vortex. We also investigate a case not covered by the experiments or the previous numerical works: 
for strong quartic potentials, the giant vortex is obtained for lower Q, before the lattice is formed. We analyze 
in detail the three dimensional structure of vortices. 

PACS numbers: 03.75.Fi,02.70.-c 
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I. INTRODUCTION 

The existence and formation of quantized vortices have re- 
cently been widely studied in Bose Einstein condensates [1- 
4, 6, 7]. One type of experiments consists in rotating the mag- 
netic trap confining the atoms. For a harmonic trapping poten- 
tial {1 / 2)mui 2 1 _r 2 , and a rotating frequency fl close to 0.7to±, 
vortices start to appear and arrange themselves into a lattice 
[5]. As ft is increased, the number of vortices increases as 
well. In the case of a harmonic trap, the confinement and 
the centrifugal force prevent the condensate from rotating at a 
frequency 57 beyond u>±. The regime of fast rotation is espe- 
cially interesting since it provides a setting for a large number 
of vortices and eventually giant vortices [8, 9]. 

Theoretical and numerical studies have considered stiffer 
potentials than the harmonic one, behaving like r™ or r 2 + r 4 
[10-13]. This type of trapping, which eliminates the singu- 
lar behavior at O = u>±, has recently been achieved experi- 
mentally by superimposing a blue detuned laser beam to the 
magnetic trap holding the atoms [14]. The resulting potential 
is 



with 



1 



V trap {r,z) = V h (r,z) + U(r), 
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Vh = -muj_r + -muo z z , and U(r) = J7 cxp( j 



ur 



(2) 

For r/w sufficiently small, the resulting potential can be ap- 
proximated by : 
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— T r + -muj z z . (3) 



The purpose of this paper is to find the stable states (vor- 
tex lattice, vortex array with hole and giant vortices) of the 
condensate with this type of trapping potential and to analyze 
their three-dimensional structure. We consider a case similar 



to the experiments and previous theoretical settings, where the 
amplitude U$ of the superimposed laser is small, so that the 
coefficient of the r 2 term is positive. But we are especially in- 
terested in the case where the laser beam has sufficiently large 
amplitude so that 



2E/o 



(4) 



This changes the sign of the harmonic part of the potential (3). 
The point is that, this case of a quartic minus harmonic poten- 
tial allows to observe giant vortices at lower angular velocities 
than previously and the structure of vortices is different. 



II. NUMERICAL APPROACH 

We consider a pure BEC of atoms confined in a trapping 
potential Vt rap , rotating along the z axis at angular velocity O. 
The equilibrium of the system corresponds to local minima of 
the Gross-Pitaevskii energy in the rotating frame 
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where g^u — ^h 2 a/m and the wave function <j> is normal- 
ized to unity j v \(j>\ 2 = 1. 

For numerical purposes, it is convenient to rescale the vari- 
ables as follows: r = x/R, u(r) = i? 3 / 2 0(x), where 
R = dj y/e and 
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^moj± / \8ttNo, 
In this scaling, the trapping potential (3) becomes 

V = {l-a)r 2 + ^kr i +f3 2 z 2 , 



, Q = Q/(eu±). (6) 
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where 
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Note that we take lo \ 



(which is the frequency of the origi- 

as a scal- 



nal harmonic potential Vh), and not uj± — a\ 
ing frequency for 51. For numerical applications, we choose 
e = 0.02, = w z /co ± = l/7,k/a = 0.25, which fit the ex- 
perimental values of Ref. [14]. In [14], a — 0.25, but we will 
take bigger values since our aim is to understand the influence 
of a when it gets bigger than 1 . 

Then, we use the dimensionless energy introduced in [15] 



E(u) = H(u) 
where H is the hamiltonian 



QL z (u), 



H{u) 
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and L z the angular momentum axis 

du 
ox 



L z (u) = i j 
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du 
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Using a hybrid Runge-Kutta-Crank-Nicolson scheme de- 
scribed in Ref. [15], we compute critical points of E(u) by 
solving the norm-preserving imaginary time propagation of 
the corresponding equation: 



du 1 
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u(V+\u\ 2 ) + fi £ u, (12) 



where \i e is the Lagrange multiplier for the constraint 
fjy \u\ 2 = 1 and with u = on dT> and. Here, D is a rectan- 
gular domain containing the condensate. A typical simulation 
uses a domain (x, y, z) e [-2, 2] x [-2, 2] x [-2.8, 2.8] with 
a refined grid of 200 x 200 x 140 nodes, which is sufficient 
to achieve grid-independence for all considered numerical ex- 
periments. 

We first compute the steady state corresponding to a non- 
rotating (51 = 0) condensate, using as initial condition u = 
yffh^, the Thomas-Fermi profile 

p TF (r) =p + (a- l)r 2 - ifcr 4 - 2 z 2 . (13) 

Depending on the choice of a, the Thomas-Fermi density 
profile can display three different shapes, as shown in figure 
1. The corresponding steady solutions obtained for 51 = 




FIG. 1: Thomas-Fermi limit p TF for different values of a. 

(which will be used as initial conditions for the subsequent 
runs with 51 > 0) are displayed in figure 2. We can distin- 
guish three cases: 



,0 3 . 



FIG. 2: Different shapes of the condensate at 51 = 0: isosurfaces of 
lowest density in the condensate for a = 0.9 (picture 1), 1.1 (picture 
2), 1.2 (picture 3). 



• a < 1 (weak quartic case) is the case closest to the ex- 
periments and is strongly influenced by the harmonic 
part. For fl = 0, a classical prolate condensate is ob- 
tained. As 51 increases, the effective trapping potential 
V e ff(r) = V(r) — e 2 51 2 r 2 starts to have a mexican 
hat structure. A vortex lattice appears for intermediate 
values of 51 and turns into a lattice with a hole for large 
51. 

• a > 1 (intermediate quartic case): the density profile 
has a depletion close to the center at 51 = but no hole. 
The criterion for this case is 



6= ^ fc5/4 >1 



(14) 



The density profile starts to have a hole for intermediate 
values of 51. 

• a > 1 and £ < 1 (strong quartic case): the density 
profile has a hole for all 51. 



III. DESCRIPTION OF THE RESULTS 

Depending on the values of a and 51, we observe different 
types of configurations: vortex free configurations where the 
amplitude of the wave function takes into account the shape of 
the effective trapping potential, vortex lattices, vortex arrays 
with hole and giant vortices. 



A. Intermediate quartic case (a — 1.1) 

The potential V has a Mexican hat structure. The isosurface 
of lowest density of the solution is plotted in figure 3, the top 
view in figure 4 and in the middle plane z — in figure 5. 
For 51 small, the density has a depletion close to the center 
of the condensate but no hole and no vortices. For 51 larger 
(fl/u>± > 0.16), vortices are nucleated. 

For 0.16 < fl/cLi± < 0.24, the density of the solution is 
zero close to the top and bottom of the condensate, but not at 
the center, which gives rise to a special structure of vortices: 
the vortices arrange themselves along two concentric circles. 
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FIG. 3: (a = 1.1) Side view of the condensate for £l/ui± — 0.12 
(a), 0.2 (b), 0.28 (c), 0.32 (d). Isosurface of lowest density. 
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FIG. 4: (a = 1.1) Top view of the condensate for fi/u>j_ = 0.12 (a), 
0.16 (b), 0.2 (c), 0.24 (d), 0.28 (e), 0.32 (f), 0.4 (g) and 0.48 (h). 
Isosurface of lowest density. 
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FIG. 5: (a = 1.1) Density contours in the plane z — for the same 
states as in figure 4: Q/uj± = 0.12, 0.16, 0.2, 0.24, 0.28, 0.32, 0.4 
and 0.48. 



The inner circle is made up of vortices which are isolated in 
the center of the condensate but reconnect towards the top of 
the condensate (see the details in figure 6). The outer circle 
is made up of almost straight U vortices that reconnect to the 
inner circle close to the top and bottom of the condensate. As 
increases, the number of vortices on each circle increases. 
In figure 4(b), the inner vortices seem to be bigger, but this is 
just an effect due to the projection and the bending: the view 
at z = (figure 5) allows to check that all vortices have the 
same size. 

For tt/uj^ > 0.24, the density profile of the solution is 
zero in the center of the condensate, hence this creates a giant 
vortex: the straight vortices that were close to the center on 
the inner circle have merged into a giant vortex. There are 
also isolated vortices regularly scattered on a circle around the 
giant vortex. As £1 increases, the number of vortices inside 
and outside the giant vortex increases and the length of the 
isolated vortices decreases as can be seen in figure 3. 
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FIG. 6: (a = 1.1) Vortex reconnection details for Q./uj± = 0.16 (a), 
0.2 (b). 

Note that the isolated vortices are U vortices that reconnect 
to the giant vortex at the center, not to the boundary of the 
condensate, as in the case of the harmonic trapping [15], that 
is their bending is concave not convex. 

For Vl I oj± = 0.48, (see figure 7) the number of vortices has 
increased and there are 2 outer circles of vortices around the 
giant vortex: one circle of U vortices that reconnect to the gi- 
ant vortex and one circle of vortices that reconnect to the outer 
boundary of the condensate. Both have different concavity in 
their bending as illustrated in figure 7. 




FIG. 7: (q = 1.1) Vortex details for n/co ± = 0.48. 



B. Strong quartic potential case (a = 1.2) 

The effective potential has a Mexican hat structure for all ft 
and the density profile of the solution always has a hole in the 
center as illustrated in figure 8. 




FIG. 8: (a — 1.2) Top and side view of the condensate for Q/luj_ — 
0(a), 0.12(b) and 0.2(c). 



For small 0, there are no vortices, that is L z = 0; it is only 
the modulus of the solution that goes to zero. For larger 
(£l/u± > 0.12), the hole contains a giant vortex and L z in- 
creases with ft (see figure 9). We have not found any isolated 
vortex around the giant vortex: all vortices are included in the 
central giant vortex because of the strong potential. The giant 




FIG. 9: Angular momentum L z (in units of h) for all studied config- 
urations. 

vortex phase profiles (figure 10) show that the phase singulari- 
ties do not completely overlap in the center of the vortex. This 
feature has already been observed in two-dimensional numer- 
ical simulation of a fast rotating condensate by Kasamatsu et 
al [13]. They described the giant vortex as the hole containing 
single quantized vortices with such low density that they are 
discernible only by the phase defects. 
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FIG. 10: (a = 1.2) Phase distribution in a central z = cut plane. 
fi/w_L =0.12 (a) and 0.2(b). 



C. Weak quartic case (a = 0.9) 

This is the case closest to the experiments [14]. The special 
feature of this case is that one has to achieve larger values of 
il in order to obtain giant vortices. The density profile of the 
solutions are shown in figures 1 1 and 12. 
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FIG. 11: (a = 0.9) Top view of the condensate (up) and density 
contours in the plane z = (down) for Q/ui± — 0.32 (a), 0.4 (b), 
0.48 (c) and 0.56 (d). 



Figure 12 show the three dimensional structure of vortices. 
There are isolated single quantized vortices, forming a lattice. 
Increasing leads to a denser lattice (20 vortices for £1 /oj± = 
0.4 and 38 for fl/ui± = 0.48). As a consequence, the angular 
momentum (figure 9) grows rapidly to high values. 

^From Q/uj± = 0.48, the vortices near the center of the 
condensate start to merge, leading to a central structure similar 
to that displayed in figure 6(a). For n/oj± > 0.56, the central 
vortices have merged into a giant vortex. The lattice still exist 
around. Similarly to the experiments, the hole is obtained for 
large values of angular velocity (fl/u;± > 0.56). 

It is interesting to note from the side view of the conden- 
sate (figure 12) that most vortices of the lattice are straight, 
but some bent vortices (U shape) exist. The U vortices are 
either connected to the outer boundary (bending outwards) of 
the condensate (figure 12a,c) or to the giant vortex (bending 
inwards) (figure 12d). 
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FIG. 12: (a = 0.9) Side view of the condensate for Q/uj± = 0.32 
(a), 0.4 (b), 0.48 (c) and 0.56 (d). 



IV. CONCLUSION 

We have studied stable configurations of the Gross 
Pitaveskii energy when the trapping potential is modified to 
include a quartic minus a harmonic term. 

For weak quartic potentials, the solution evolves from a vor- 
tex lattice to a vortex array with hole when the angular veloc- 
ity f2 is increased. For stronger quartic potentials, giant vor- 
tices are obtained for lower O, at a stage where the lattice is 
not so dense. The typical structure of vortices is to have a cen- 
tral giant vortex with an outer circle of vortices around. We 
believe that there should be a criterion depending on the radius 
of the condensate and the radius of the annulus that should 
characterize the final structure of the giant vortex: whether 
there is or not a circle of vortices around the giant vortex and 
its precise location. 

The form of the potential considered in our simulations was 
inspired from recent experiments [14]. We have checked that 
keeping the exponential part of the potential instead of its 
quartic minus harmonic approximation does not change the 
qualitative behaviour of the solutions. This suggests that if 
this situation could be achieved experimentally, it would allow 
to observe giant vortices for lower velocities than previously, 
that is before reaching the fast rotation regime. 
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